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Abstract 

A simple model allows us to study the nonclassical behavior of slowly moving atoms 
interacting with a quantized field. Atom and field become entangled and their joint 
state can be identified as a mesoscopic "Schrodinger-cat" . By introducing appropri- 
ate observables for atom and field and by analyzing correlations between them based 
on a Bell-type inequality we can show the corresponding nonclassical behavior. 
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1 Introduction 

In 1935 Erwin Schrodinger published an influential paper [1] whose contents 
is still a timely topic. In this work he introduced the term entanglement which 
lies at the very heart of quantum theory but has no counterpart in the clas- 
sical description of nature. The measurement result found for one part of an 
entangled system depends in a nonclassical way on what is measured on the 
other part. This — in a certain sense — nonlocal feature of quantum me- 
chanics has been pointed out first by Albert Einstein, Boris Podolsky and 
Nathan Rosen (EPR) in another famous paper of the year 1935 [2]. In fact, 
Schrodinger himself called his paper a "confession" sparked by the EPR work. 
It was not until 1964 when John Bell succeeded to formulate the consequences 
of the EPR propositions on locality and reality [3]. The limitations for specific 
correlations imposed by a classical theory based on locality and reality found 
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expression in his seminal Bell-inequality. Their violation therefore provides a 
clear signature of nonclassical behavior of entangled systems [4]. 

The incompatibility between quantum mechanics and a classical description 
becomes also evident when one extrapolates the quantum mechanical descrip- 
tion to macroscopic objects. Schrodinger highlighted this fact also in his paper 
of the year 1935 when he introduced his famous gedankenexperiment of an iso- 
lated cat in a box [1]. The formulation of his gedankenexperiment results in 
the question whether superpositions like 



exist in nature or if they are just a mathematical consequence of quantum 
mechanics. Here, the states \ alive) and \dead) mean macroscopic states of 
the "Schrodinger cat" correlated with the exited state | e ) and the ground | g ) 
of an atom [1]. Since the work of Schrodinger it has become clear that such 
states can be at least prepared in the laboratory [5,6,7,8,9,10] and possess the 
properties predicted by quantum theory. 

In this paper we present an atom-optic model which allows us to create at 
least mesoscopic Schrodinger-cat states. The cat-like degree of freedom will be 
replaced by the motion of an atom. Likewise we will point out the nonclassical 
aspect of this state by analyzing a Bell-type correlation sum. Hence it is the 
aim of the present paper to combine [11,12,13] Schrodinger's cat with his term 
entanglement using only basic atom-optic elements [14,15]. The used model 
is simple and only applies operationally defined elements. 

The paper is organized as follows. In Sec. 2 we introduce a simple model 
consisting of a single two-level atom moving slowly through a quantized field. 
We show that the time evolution of the combined system leads to a state very 
similar to the original "Schrodinger-cat" state. In Sec. 3 we show how one can 
construct observables that demonstrate nonclassical correlations between the 
quantized atomic motion and the cavity field. Finally, in Sec. 4 we address 
the question if and in which way the observables constructed in the previous 
section have an operational meaning. We conclude with Sec. 5. 



2 Atomic Schrodinger cat 

The aim of the following section is to briefly review a simple model based 
on atom optics which allows us to create states very similar to the original 
Schrodinger-cat state, Eq. (1). 




(1) 
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We consider [16,17] a cold two-level Rydberg atom with ground state | g ) and 
excited state | e ) moving slowly in z-direction through a single-mode quan- 
tized field inside a cavity, as shown in Fig. 1. The atom interacts resonantly 
with the cavity field which is initially cooled down to the vacuum state | 0). 
We describe the initial center-of-mass motion of the atom by a Gaussian wave 
packet 

(z\&) = $(z) = (27ra)- 1/4 exp 



Aa 2 



+ ikoz 



(2) 



of width a, centered at — z . At time t — the wave packet starts moving with 
mean velocity hko towards the cavity "from the left to the right" . 

Let us assume that atom and field are initially uncoupled. When we start with 
an excited atom, the initial state of the complete system reads 

|tf(t = 0)) = |$)<g>|e)<g>|0). (3) 



In the interaction picture the time evolution of the combined system is de- 
scribed by the Hamiltonian 

H=p- + H'yu(z)(a-a i + a + a). (4) 
2m v ' 



The atom of mass m enters the cavity mode whose spatial structure is de- 
scribed by the mode function u(z). The resonant interaction inside the cavity- 
field is modeled by a typical Jaynes-Cummings term [20] with coupling strength 
hq. The operator-combination a_a) describes the decay of an atom and the 
simultaneous creation of a cavity mode photon and has to be superposed with 
the combination a + a standing for the excitation of the atom and the annihila- 
tion of a cavity mode photon. It turns out that exactly this resonant exchange 
of excitations between two-level atom and cavity field leads to a spatial split- 
ting of the initial atomic wave packet | $ ) . Moreover, it leads to nonclassical 
entanglement between atom and field. 

To see this we look at the time evolution 

| V(t) ) = exp 

of our system in terms of the so-called dressed states 

\n,±) = ^={\e)\n-l)±\g)\n)}, 



$)|e)|0) 
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which combine the internal states of the atom with the Fock states | n ) of 
the cavity field. Since the dressed states are the eigenstates of the Jaynes- 
Cummings part of the Hamiltonian, Eq. (4), with eigenvalues ±\/n, we can 
appropriately expand the initial state \^/(t)), Eq. (3) and arrive at 

= ^[\& +] (t))\l, + ) + \&-\t))\h-)}, (5) 



where we have introduced the partial wave packets 



$ [±1 (t)) = exp 



H \ 2m 



± h'yu(z) t 



Hence the complete time evolution of \^(t)) can be reduced to the time evo- 
lution of the partial wave packets | ) in the effective optical potentials 
±fryw(;z). Note that the atom "sees" both potentials at the same time. 

For a given mode function u(z) the partial wave packets | $ [±] (t) ) can now be 
calculated. We concentrate here on the most simple case, that is we assume a 
square-well u(z) = \z\) of length 2L for the envelope of the cavity-field. 

It has been shown in [16] that despite an additional phase the Gaussian form 
of the initial wave packet persists for both partial wave packets. On the right 
hand side of the cavity, that is for times t > T, where 



T = 



2z n 



hk /m 



represents the characteristic time needed by a classical particle moving form 
— z to +z with constant momentum hk , the partial wave packets read 



(z\$W(t)) = $ [±] (z,t) = 



e i[Tk D/2+<p(z,t;z ±D/2)] 



exp 



(z + z - ± D/2) 2 



i(a 2 + 



h 2 t 2 

Aa 2 m 2 , 



.(6) 



Here we have introduced the phase 



ip(z, t; z ) = k [z- —t + 



^oA {z + z,-^tf 



2m J 4a 2 ( 



2ma 2 i ht 



ht 2mo 



The maxima of the two wave packets &^(z,t), Eq. (6), are spatially well 
separated by the distance 



V 2m ) 



D = 77^2L, (7) 
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which scales with the length 2L of the interaction zone. Note that we have 
started initially with an empty cavity: even the vacuum state strongly influ- 
ences the atomic motion and splits it into two spatially well separated parts 
[21,22]. 

When we now compare the evolved state, Eq. (5), with the original Schrodinger- 
cat state, Eq. (1), we clearly recognize the same structure. In both cases we 
have a superposition of a mesoscopically or even macroscopically distinguish- 
able state entangled with a microscopic state. The only difference is that in 
contrast to the original cat state, where cat and atom are two completely sep- 
arate objects, here the dressed states involve both the internal atomic states 
as well as the states of the cavity field. In the following we will show the non- 
local and therefore nonclassical properties of this cat-like state by analyzing 
a Bell-type sum of correlation functions. 



3 Nonlocal field-atom correlations 



The time evolution of our simple model-system presented in the previous sec- 
tion entangles all three degrees of freedom, i.e. the center-of-mass motion of 
the atom, its internal states and the cavity-field mode. Additionally, after 
passing the cavity the atom evolves freely in z-direction and is therefore in 
principle spatially well separated from the cavity. Nevertheless, measurement 
results found for the center-of-mass motion may depend nonlocally on what 
is measured on the internal states and the cavity-field mode. Such nonclas- 
sical correlations between subsystems can be shown using Bell-type inequali- 
ties [3,4,23]. Such inequalities always involve correlation functions of the form 
(A(a)B(P)) = \A(a)B(f3)\ ^ ) with hermitean and dichotomic operators 
A(a) and B((3) describing measurements on a state | ^ ) depending on the 
measurement parameters a and (3 respectively. In particular we will show that 
a CHSH-sum of correlations [23] 



B QM (d,a'JJ') = (A(d)B0)) + (A(a)B0)) 

+(A(d')B0)) - (A(d')B0)) (8) 

can lead to a regime | Bqu | > 2 for our atom-field system and therefore clearly 
reveal its nonclassical behavior. For corresponding experiments on spin-like 
systems we direct attention to Refs. [24,25,26,27] and references therein. 
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3. 1 Observables and correlation function 



On the basis of the considerations of the previous section we are now in the 
position to show the nonlocal character of our model-system. For our purpose 
it is appropriate to rewrite the time-evolved state Eq. (5) in terms of the 
original internal and field states which yields 



*(t)) = \ [(l^ [+1 W> + l^ H (0>)|e>|0) 

+ (|$W(t))-|$H(t)))|^)|i)]. ( 9 ) 



This state does not describe a discrete spin system but the continuous mo- 
tion of an atom coupled to a quantized field. It has been pointed out earlier 
[11,12,13,28,29] that continuous systems can violate Bell-type inequalities as 
well. Therefore, we face now the problem of finding proper dichotomic observ- 
ables for the atomic motion, the internal atomic states and the cavity field 
which allow us to obtain nonclassical CHSH-correlations, Eq. (8). 

One of the simplest dichotomic operators for the atomic motion, i.e. for the 
partial wave-packets | ) is the parity operator 



P = 



= J dz\-z)(z\. (10) 



As we know from Eq. (8) we need adjustable parameters a, (3, etc. in order 
to define the reference frame of the measurement. It is probably most natural 
[11,12,13] when we look at parity from different points (z,k) in phase-space 
by displacing the parity with the unitary displacement operator 



D(z, k) = exp i(kz — zk) 



where z is the position operator and k = p/h represents the scaled momentum 
operator. 

Although our idea is to show the nonclassical character of the CHSH-sum, 
Eq. (8), with the help of continuous degrees of freedom, i.e with the atomic 
center-of-mass motion and the cavity field we have to deal with the internal 
states. A simple projection to the ground or the excited state does not come 
into question since it would clearly destroy the entanglement between atomic 
motion and cavity field. Instead we consider a superposition measurement 
described by the hermitean operator 

a=\g)(e\ + \e)(g\. (12) 
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Merging the operators Eqs. (10), (11) and (12) we arrive finally at the observ- 
able 

A(z, k) = D(z, k)PD^(z, k) ® a (13) 
for the atom. 

Due to the chosen initial conditions, Eq. (3), the whole setup consisting of 
atom and field contains only a single excitation. Therefore, for the field we 
can construct as well a dichotomic hermitean operator of the form 

B(P) =e i/3 |0)(l| + e- t/3 |l)(0| (14) 
depending on the parameter (5. 

Having introduced all necessary ingredients we are now in the position to 
combine Eqs. (9), (13) and (14) and calculate the correlation function 

(A(z, k)B(P)) = ( (f) \A(z, k)B((3)\ *(*)), 

which depends on three real parameters: The pair (z, k) parameterizes the 
reference frame of our first subsystem, namely the atom, while the phase f3 
parameterizes the measurements on the second subsystem which is the cavity 
field. 

We can further analyze this correlation function 

(A(z, k)B(P)) = X - (( \W(z, k)\ ) - ( $ H \W(z, k) | $H )) cos p 

_!(($[+] | vv(z, k)\ $ H ) - ( $ H \W(z, k) | $ W )) sin ^15) 

when we introduce the notation 

W(z, k) = D(z, k)Pb\z, k) (16) 

for the shifted atomic parity according to Eq. (13). 

The computation of the correlation function, Eq. (15), is reduced to that of 
the different combinations (<J>M |VV(z, k) \ ). It is a straight-forward cal- 
culation [18,19] to show that for an arbitrary state | $ ) the operator W(z, k), 
Eq. (16) leads to an expectation value of the form 

oo 

($\W(z,k)\<S>) = J dye~ 2iky <5>\z-y)<$>(z + y) 

— oo 
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which is in fact the Wigner function [20] up to a normalization factor. A closer 
investigation for the partial wave packets ^^(z,t), Eq. (6), shows that the 
resulting Wigner functions are as well Gaussian distributions, now both in 
position- and in momentum-space. The explicit expressions appearing in the 
correlation function, Eq. (15), read 



\W(z, k)\^ l+] )=W + (z, k) 

■ exp | — 2a 2 (k 



- k f - 



2a 2 



hk 

z + z 1 

m 



D 
~2 



(17) 

x 2" 



($ H \W(z, k)\<^ l ' ] )=W(z, k) 



■ exp 



-2a 2 (k-k ) 2 - 



and for the interference term we find 



2a 2 



hk 

z + z 1 

m 




% - (( $ [+] \VV(z, k)\ $H ) - ( $H |yv(z, k)\ $ [ 



= exp 



-2a 2 (k - k f 




W mt (z, k) 



sin[D(ife-2A;o)]. (19) 



These functions remain stationary around the momentum k = k while in 
position space they move with velocity {hk)/m in positive z-direction and 
additionally get sheared. Moreover, each of the two Wigner functions W + (z, k) 
and W~(z, k), Eqs. (17) and (18), is shifted by half of the distance D, Eq. (7), 
originating from the interaction with the cavity field. In contrast to this the 
interference term W mt (z, k), Eq. (19), remains centered but its momentum 
envelope oscillates with frequency D. 

The correlation function (A(z,k)B(/3)), Eq. (15), therefore has a clear inter- 
pretation in phase space: Fig. 2 shows the typical picture of a Schrodinger-cat 
in phase space which results from Eq. (15) for equal-weighted trigonometric 
functions, i.e. for f3 — 7r/4. However, due to the special structure of the corre- 
lation function here the wave packet W~ (z, k) appears negative, whereas the 
wave packet W + (z, k) is always positive. 

Additionally to the measurement parameters z, k and f3 the momentum hk , 
the position — zq and the width a of the initial atomic wave packet, Eq. (2), 
enter in the quantity (A(z,k)B(/3)). Moreover, the measurement time t > 
T, the coupling strength 7 between atom and field and the length of the 
interaction zone 2L can be varied in a certain range. 
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By scaling all lengths with half the width L of the interaction zone 2L, all 
wave numbers by yj 2m^ / K and additionally fixing the product \J 2m^ / h ■ L of 
the two scaling-parameters, the correlation function, Eq. (15), depends in all 
on eight parameters. 



3.2 Nonclassical behavior of the combined atom-field system 



The final step in order to identify the nonclassical regime of the CHSH-sum 
of correlations, Eq. (8), is now to combine four of these correlation functions, 
Eq. (15), that is 

B QM = (A{z,k)B{(3)) + (A(z,k)B((3')) + (A(z',k')B((3)) - (A(z', k')B((3')).(20) 



The CHSH-combination Bqm depends on six real parameters: the positions 
z and z' , the wave numbers k and k' and the phases f3 and f3'. These six 
parameters have to be varied in order to find nonclassical values Bqm > 2, 
whereas the initial values of the wave packet, Eq. (2), and the time t > T of 
the measurement are fixed for all four correlations. 

In our case the CHSH-combination, Eq. (20) can be written as 

Bqm =X 1 -cosP + Y 1 - sin (3 + X 2 ■ cos/3' + Y 2 ■ sin (3' (21) 



where the Xi and Yi depend on z,z',k and k'. 

This form allows us to maximize Bqm with respect to the phases (5 and /3'. 
The full expression of Bq M = max£>QM can then be found by inserting the 

explicit expressions of the wave packets, Eqs. (17) - (19). We find 



B' 



QM 



-2d' 



(e " 2 x1 sinh[2a;<i] + e K ' 2 x ' 2 sinh[2a;'<i 



+ (e-* 2 -* 2 sin[2rf(/t - «(,)] + e -' 2 -' 2 s in[2rf(/ ^ , - kq)])' 



+ 



e- 2d2 (e-* 2 -* 2 sinh[2xrf] - e - K ' 2 - x ' 2 sinh[2x'rf])' 



+ (eT* 2 -* 2 sin[2d(« - ko)] - e ^' 2 -' 2 sin[2d(«' - Kq)])' 



where we have introduced the abbreviations x — (z + Zq — hkt/m) / (y/2a), k = 
\/2a{k — ko), d = D / (2y/2a) and k = \^2ak and the analogous expressions 
for the primed variables. 
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A numerical analysis of this expression shows the following result: The max- 
imal value of Bqm can be found for k — k' — and for x = —x' = 0.371. In 
other words, the nonclassicality is maximal for wave numbers k = k' = k , i.e. 
for the initial momenta Hk , and for positions x and x' displaced symmetrically 
around zero. Inserting this values and optimizing the remaining parameters 
tod= 0.741 and « = (2n+ 1) • 7r/(4d) = (2n + 1)tt/2.954, we find a maximal 
value B^gf = 2.324. 

Fig. 3 shows the CHSH-sum as a function of the parameters z'/L and f3' for 
k — k' — ko, while we optimized the values of z/L and f3 numerically. The 
region with values Bqm > 2 is dyed and shows for the chosen parameters a 
clear nonclassical regime of the CHSH-sum with a maximum of Bqm ~ 2.324. 

The numerical optimization of Bqm shows that the maximal value 2v^2 allowed 
by quantum theory [32] cannot be achieved. It seems that the special choice of 
the atomic observable A(z,k), Eq. (13), leading to Gaussian-shaped correla- 
tion functions, is responsible for this fact: They smooth the CHSH-correlation, 
Eq. (20), and therefore reduce the maximal possible value of 2\/2. 



4 Operational meaning of the observables 

In the previous section we have introduced the operators A(z,k) and B(f3), 
Eqs. (13) and (14), for the atom and the field. A combination of four correlation 
functions has led us to the CHSH-correlation sum, Eq. (20), which was shown 
to reach nonclassical values for a specific choice of the parameters of our 
system. This implies the possibility of a simultaneous measurement of the 
atomic operator A and the field operator B. In the present section we shall 
shortly present the basic ideas how each of the occurring operators can be 
realized experimentally. 

We concentrate first on the atomic operator A(z, k), Eq. (13), and recall that 
it consists of two parts, namely of the operator W(z,k), Eq. (16), for the 
atomic center of mass motion and the operator a, Eq. (12), for the internal 
atomic states. 

The suitable operator VV(z, k) for the atomic motion turned out to be the par- 
ity, displaced in phase-space by the amount (z,k). Thereby, the phase-space 
point (z, k) had to be chosen such that the CHSH-correlation sum moved to 
the nonclassical regime. In order to realize this operator experimentally we 
first recall how the expectation value ( | VV(0, 0) | $^ ) of the parity oper- 
ator could be measured. Since consecutive energy eigenstates of a symmetric 
potential have alternating parity we could for example determine the parity 
of the atomic wave packet by bringing it into a symmetric potential and mea- 
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suring its energy. However, we are not interested in the parity of the atomic 
wave packet at the origin but at the phase-space point (z, k), determined by 
the displacement operator D(z, k), Eq. (11). This operator acts on the atomic 
wave packets | ). It is equivalent to an instantaneous displacement of the 
atom from the origin to the point (z, k). We can as well imagine to displace 
the potential by the same amount in the opposite direction, which might be 
experimentally easier than displacing the atom. 

The operator a, Eq. (12), for the internal atomic states has the structure of 
a Pauli-spin matrix and describes a measurement of the atomic dipole. The 
experimental realization can be achieved by state-selective field ionization of 
the atoms [30]. 

For the field measurement we take into account that we have at most one exci- 
tation in the cavity. Therefore, the operator B{(3), Eq. (14), also has the struc- 
ture of a Pauli-spin matrix, but with an additional phase. The experimental 
realization of this operator is equivalent to an electric-field-measurement of 
the cavity. Such measurements have been done by state-of-the-art homodyne 
techniques (see e.g. [31]). In principle, the field is coupled out from the cavity 
and due to a beam splitter becomes entangled with a coherent state | |a|e l/3 ) 
of large amplitude. The differences of the intensities at the two outcomes of 
the homodyne interferometer finally determine the electric field of the cavity. 
The phase (5 is thereby determined by the phase of the coherent state. 



5 Conclusions 



In this paper we have proposed a simple dynamical model and presented how 
its entangled continuous degrees of freedom could be used to exhibit the nonlo- 
cal aspect of quantum mechanics. We have investigated the operators suitable 
in this context and have given ideas how to define them operationally. In a 
further publication, we were able to show [33] that for atoms in a linear po- 
tential, e.g. under the influence of gravity, this nonlocal feature remains still 
valid. However, there might exist other, maybe simpler operators which would 
indicate nonlocality as well. Moreover, it is not yet clear, if other criteria for 
nonclassical correlations [34,35,36] can be applied to our model. The influence 
of decoherence effects [37] is another topic we have not yet taken into account 
in this context. 

This work was financially supported by the Deutsche Forschungsgemeinschaft 
and by the European Community's Human Potential Programme under con- 
tract "QUEST" HPRN-CT-2000-00121. K.W.'s work was additionally par- 
tially supported by a KBN grant 2P03 B 02123. 
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Fig. 1. The Gaussian wave packet of a two-level atom in the excited state moves 
freely towards a cavity. The field inside the cavity is cooled down in the vacuum 
state. The atomic two-level structure and the quantized field create two effective 
optical potentials which split the wave packet in a superposition state. 




Fig. 2. Correlation function Eq. (15) in phase-space. For a equal weight of the 
two Gaussian distributions W ± (z, k) and the interference term W mt (z, k), i.e. for 
/3 = 7r/4, we get the typical picture of a Schrodinger-cat. For the generation of this 
plot we scaled all length with half the width L of the resonator field, all wave num- 
bers with (2my / K) 1//2 and the time with 7, where m is the mass of the atom and 7 the 
coupling-strength between atom and cavity field. The plot results from the param- 
eter settings (2m 7 //i) 1 / 2 L = 250, z /L = 1.75, a/L = 0.03, k /(2m- f /h) 1 / 2 = 2.31. 
and 7^ = 190. The spatial separation between W + (z, k) and W _ (z, k) is therefore 
given by D/L = 2(2m-f/h)/k% = 0.375. 
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Fig. 3. Nonclassical region of the CHSH-correlation sum, Eq. (20). The CHSH-sum 
£>qm is plotted over z'/L and (3' while z/L = 0.24 and (3 = ir/2 are fixed. The param- 
eters are (2mj/H)^ 2 L = 6030.7, z /L = 0.15, a/L = 0.067, k /(2mj/h) 1 ^ = 2.668 
and jt = 400. The dyed region shows the nonclassical regime. For this parameters 
the maximal value is given by Bq^ ~ 2.324. 
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